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Abstract. For the root system of each complex semi-simple Lie algebra of rank two, 
and for the associated 2D Toda chain £ = {u xy = exp(Ku)} , we calculate the two 
first integrals of the characteristic equation D y (w) = on £. Using the integrals, we 
reconstruct and make coordinate-independent the (2 x 2)-matrix operators □ in total 
derivatives that factor symmetries of the chains. Writing other factorizations that 
involve the operators □, we obtain pairs of compatible Hamiltonian operators that 
produce KdV-type hierarchies of symmetries for £. Having thus reduced the problem 
to the Hamiltonian case, we calculate the Lie-type brackets, transferred from the 
commutators of the symmetries in the images of the operators □ onto their domains. 
With all this, we describe the generators and derive all the commutation relations in 
the symmetry algebras of the 2D Toda chains, which serve here as an illustration for 
a much more general algebraic and geometric set-up. 



Introduction. In the paper [8], we introduced a well-defined notion of linear matrix 
operators in total derivatives, whose images in the Lie algebras of evolutionary vector 
fields on the jet spaces are closed with respect to the commutation. This yields a 
generalization for the classical theory of recursion operators and Poisson structures for 
integrable systems. We explained how each operator transfers the commutation of the 
vector fields to the Lie brackets with bi-differential structural constants on the quotient 
of its domain by the kernel. 

Second, we associated such operators with the 2D Toda chains 

m 

^Toda = [ui y = exp(^ K)u j ) , 1 < i < m} (1) 

3=1 

related to semi-simple complex Lie algebras [TU1 fTT] . We derived an explicit formula 
for the operators that factor higher symmetries of these chains. Using the auxiliary 
matrix operators that we proved to be Hamiltonian, we elaborated a procedure that 
yields all the commutation relations in the symmetry Lie algebras sym ^Toda (naturally, 
these symmetry algebras are not commutative). This solved a long-standing problem in 
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geometry of differential equations and completed previously known results by Leznov, 
Meshkov, Shabat, Sokolov, and others (see [HI Q21 [ESI 20J and references therein). 

Actually, the general scheme of [8] is applicable, in particular, for the description of 
symmetry algebras for a wider class of the Euler-Lagrange hyperbolic systems of Liou- 
ville type [161 ED] • Moreover, the group analysis of integrable systems, as a motivation, 
results in the well-defined concept of operators whose images span involutive distribu- 
tions on the jet spaces, but not on differential equations, which is of an independent 
interest. 

In this note, we illustrate the reasonings of [8 J using the root systems of the complex 
semi-simple Lie algebras of rank two. Among all two-component exponential nonlinear 
systems (pQ), these 2D Toda chains with K Cartan matrices admit the largest groups of 
conservation laws [TO] and are integrable in quadratures [TO] . 

The equality of the rank to two means the following: 

• The hyperbolic Toda chains (CEJ) upon u 1 , u 2 are, we repeat, two-component. 

• The number of vector fields Yi that generate the characteristic Lie algebra 
through commutators (see section [2~Tl below and [TT]) equals two. Also, the num- 
bers of linear independent iterated commutators K^...^) = . . [Y^^YiJ . . .]] 
fall at most twice, and the accumulated sum of these dimension's falls equals 
two. Thence, by the Frobenius theorem, two invariants w 1 , w 2 appear [] Using 
the characteristic Lie algebras, we introduce two finite sequences of the adapted 
coordinates, which simplifies the description of these invariants. On the other 
hand, we use the two invariants for replacing the derivatives of the two depen- 
dent variables at all sufficiently high differential orders. 

• Conservation laws for Toda system ([!]) are differentially generated (up to x <-> 
y) by these two invariants, which are solutions of the characteristic equation 
D y (w) = on £ Toda . 

• Higher symmetries of the Toda chain ([T]) have a functional freedom and are 
parameterized by two functions 1 , <fi 2 that depend on x and any derivatives of 
the integrals w l up to a certain differential order. 

• The differential operators □ that yield symmetries of (pQ), when applied to the 
tuples (d) 1 ,^ 2 ), are (2 x 2)-matrices. 

• The Lie algebra structures transferred from sym£<roda to the domains of □ are 
described by the bi-differential brackets {{ , }}□ that contain two components. 

• The Hamiltonian structures that are defined on the domains of the oper- 
ators □ but take values elsewhere (in the Lie algebra of velocities of the inte- 
grals w\ see [S]), are also (2 x 2)-matrices. Likewise, the brackets {{ , }}A k trans- 
ferred onto the domains of from the commutators of Hamiltonian vector fields 
in their images are also two-component (thence the equality {{ , }}i fc = {{ , }}□ 
makes sense). 

• The KdV-type hierarchies of velocities of w % and the modified KdV-type hier- 
archies of commuting Noether symmetries of ^Toda, which are related by two- 
component Miura's substitutions with the former, are composed by the (right- 
hand sides of) two-component evolutionary systems. 

^For example, the paper |19j contains a brute force classification of integrable one-component hy- 
perbolic equations with respect to the low-dimensional characteristic Lie algebras. 
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All constructions and notation follow [E] except for the characteristic Lie algebras 
that were introduced in [TT] and were discussed in detail in [15J. All notions from the 
geometry of PDE are standard (see [21 El EE]) and have been surveyed in [HI sect. 2]. 
All extensive calculations were performed using the software [12] . 

To start with, we recall that in the fundamental paper [16j, A. B. Shabat et al. proved 
the existence of maximal (r = f = m) sets (j2J) of conserved densities 

w±, . . . , w r e ker D v \ c , wr, ker D x \ c . (2) 

y I tToda I ^Toda V ' 

for ^Toda if and only if the matrix K in (OQ) is the Cartan matrix of a root system for 
a semi-simple complex Lie algebra of rank r, which is always the case in what follows 
with r = 2. Note that the integrals §2§ allow to replace the derivatives of unknown 
functions of any sufficiently high order using the derivatives of the integrals. In [T5] . 
A. B. Shabat proposed an iterative procedure that specifies an adapted system of the 
remaining lower-order coordinates and that makes linear the coefficients of the linear 
first-order characteristic equation D y {w) = on £xoda- That algorithm is self-starting, 
simplifies considerably the search for the first integrals of the characteristic equation, 
and gives the estimate for the differential orders of solutions. 

Another method (which we do not use here) for finding the first integrals is based on 
the use of Laplace's invariants, see [20]- The authors of that paper investigated (pri- 
marily, in the case of one unknown function and one equation £ upon it) the operators 
that factor symmetries of £. Also there, the pioneering idea to study the operators 
whose images are closed under the commutation was proposed. We indicate further the 
papers [31 [3, [TF1 [18] that address the problem of construction of such factoring operators 
for multi-component hyperbolic systems of the Liouville type. 

The general concept of operators whose images determine involutive distributions, the 
definition itself and the classification, has been elaborated in [8]. There, as a by-product, 
we obtained an explicit formula for the operators □ that factor higher symmetries 
of the Euler-Lagrange Liouville-type systems and for the bi-differential brackets on 
their domains. The former yields all the generators of the symmetry algebras for such 
systems, and the latter describes all the commutation relations. 

This paper is structured as follows. First we outline the basic concept using the 
scalar Liouville equation u xy = exp(2w) as the motivating example. This covers the 
case of the root system Ai. The following fact, which holds true for any rank r > 1, 
is very convenient in practice: the differential orders of the r integrals w 1 , . . ., w r for 
the 2D Toda systems (TjQ) associated with the complex semi-simple Lie algebras g are 
equal (up to a shift by +1) to the gradations for the principal realizations of the basic 
representations of the respective untwisted affine Lie algebras (see, e.g., the list 
in §14.2]). 

Then we realize the geometric scheme 

D y (w) = .— > w\w 2 .— > □ .— ► A k .— > {{,}}□; A k .— > A[\ 

for the simple complex rank two Lie algebras (for the root systems A 2 , B 2 — C2, and G 2 ). 
In other words, we associate the operators to invariants, the brackets to operators, and 
find the deformations of the Poisson structures. Only once, for the root system A 2 , we 
calculate the characteristic Lie algebra for the corresponding 2D Toda chain (CQ) and 
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obtain the integrals w 1 ^ 2 using the adapted system of coordinates. We modify the 
scheme of [IS] such that, first, the two-component Toda chain is not represented as a 
reduction of the infinite chain and, second, we do not introduce an excessive third field 
which is compensated by a constraint (as in [15]). 

Remark 1 . We do not of course re-derive the structures for the algebra D2 = Ai © Ai 
that is not simple, since the operator □ and KdV's second Hamiltonian structure 
are known for each of the two uncoupled components of the chain (CQ) with K — 
see Example [1] below. However, the "x-component" of the full symmetry algebra with 
the generators 

9 V U* 2 ) \^{x,[w# l ],[w# 2 })) 

is not just the direct sum of the two symmetry subalgebras for the two Liouville equa- 
tions. Indeed, this formula shows that the integrals can be intertwined in the generators, 
although the fields are not coupled in system (CQ) with the choice of K as above, and it 
proves that the symmetries intertwine the fields. 

Remark 2. Each of the operators □, which we obtain from the r integrals, consists 
of r columns, one column for each integral. In the r = 2 case, only the respective first 
columns were specified in the encyclopaedia [I], see also [13]. Here, we complete the 
description of the symmetry generators. 

1. Basic concept 

Let us begin with a motivating example. 

Example 1 (The Liouville equation). Consider the scalar Liouville equation 

£uou = {u xy = exp(2w)}. (3) 

The differential generators w, w of its conservation laws [rj] = J f(x, [w]) dx+ J f(y, [w]) dy 
are 

w = u 2 x - u xx and w = u 2 y - u yy (4) 

such that D y {w) = and D x {w) = by virtue (=) of ^Lkm and its differential conse- 
quences. The operators 

□ = u x + \D X and □ = u y + \D y (5) 
factor higher and Noether's symmetries 

„ = □««,[«.])), ^ = n(^M) ; ?=0(*»,m», ?£ = □(«!>) 

of the Euler-Lagrange equation © for any smooth 0, and H, H. Note that the 
operator □ = \D~ 1 o Uw^Y is obtained using the adjoint linearization of w, and 
similarly for □. 

Each of the images of © is closed w.r.t. the commutation such that 

[□(?),□(<?)] =□([?, </]□), 
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where the bracket [, ]□ on the domain of □ admits the standard decomposition (the 
vector field 6^ = J2k -^(v 3 ) ' 9/duk is the evolutionary derivation), 

[p, q]n = £□(?)(?) - £□(<?) (p) + {{p, <?}}□• 

For the operator □ on the Liouville equation, the bi-differential bracket {{ , }}□ is 

{{P, q}}n = D x (p) ■ q - p ■ D x (q), 

and similar formulas hold for the operator □. The symmetry algebra sym^Lkm — 
im □ + im □ is the sum of images of (jSJ), and the two summands commute between each 
other, [im □, im □] =0 on ^Liou- Therefore, 

[im □ + im □, im □ + im □] C im □ + im □. (6) 

The operator □ factors higher symmetries of the potential modified KdV equation 

£pmKdV = {Ut = ~\u xxx + U% = □(»}, (7) 

whose commutative hierarchy is composed by Noether's symmetries (pc G im(Do 5 /5w) 
of the Liouville equation @. The operator □ factors the second Hamiltonian structure 
B2 = □ o A\ o □* for £ pm Kdv, here A\ = D^ 1 = A^ 1 is the first Hamiltonian operator 
for the potential KdV equation and equals the inverse of the first Hamiltonian operator 
for KdV. 

The generator w of conservation laws for f^iou provides the Miura substitution (jlj) 
from £ pm Kdv t° the Korteweg-de Vries equation 

^Kdv = {w t = -\w xxx + 3ww x }. (8) 

The second Hamiltonian structure for ^Kdv is factored to the product A 2 = D'o^oD, 
where B\ = D x is the first Hamiltonian structure for the modified KdV. The bracket 
{{ j }}□ on the domain of □ is equal to the bracket {{ , }}^ induced on the domain of the 
operator A2 (which is Hamiltonian and hence its image is closed under commutation) 
for <?Kdv- In what follows, we refer to these correlations as standard, see [8]. 

Definition 1 ([20]). A Liouville -type systen^S^ is a system {u xy = F(u,u x ,u y ; x,y)} 
of hyperbolic equations which admits nontrivial first integrals 

W\, w r G ker D y \ £ ^; Wi, w f G ker D x \ £ ^ 

for the linear first order characteristic equations D y \ £ (wi) = and D x \ £ (iBj) = that 
hold by virtue (=) of Si,, and such that all conservation laws for Si are of the form 
/ f(x, H) dx® J g(y, [w]) dy. 

Example 2. The m-component 2D Toda chains (0Q) associated with semi-simple com- 
plex Lie algebras [10J constitute an important class of Liouville-type systems, here 
u = (it 1 , . . . ,u m ). Further on, we consider these exactly solvable systems, bearing in 
mind that the reasonings remain applicable to a wider class of the Euler-Lagrange 
Liouville-type systems Si. 



'There exist other, non-equivalent definitions of the Liouville type systems. 
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Remark 3. The 2D Toda systems ([T]) are Euler-Lagrange, with the Lagrangian density 
L = —^(ku x , u y ) — Hl(m; x, y). The (m x m)-matrix k with the entries 

2{a i ,a j ) J_ . 

™ij I 12 1 12 I |2 j 

is determined by the simple roots of the semi-simple Lie algebra. 

Let m = dL/du y be the momenta, then it can be readily seen that the integrals to 1 , 
. . . , w m of the characteristic equation are differential functions w l = w l [m] in m. 

Proposition 1. The differential orders of the integrals w l with respect to the mo- 
menta m for the 2D Toda chains (CQ) associated with complex semi-simple Lie algebras g 
coincide with the gradations for the principal realizations of the basic (i.e., simplest non- 
trivial highest weight, see [6]) representations of the corresponding untwisted affine Lie 
algebras 

The integrals w % for a nonlinear Liouville-type hyperbolic system can be obtained 
using an iterative procedure that is illustrated in section 12.11 below. In the meantime, 
we assume that the integrals are already known. Let them be minimal, meaning that 
/ G ker D y \ £ ^ implies / = f(x, [w]). 

Theorem (|8J). Let the above assumptions and notation hold. Introduce the operator 

□=(W. (9) 
which is the operator adjoint to the linearization (the Frechet derivative) of the inte- 
grals w w.r.t. the momenta m. Then we claim the following: 

(i) All (up to x y) Noether symmetries ipc of the Lagrangial £ for £l are 

ip c = □ ^ ^ J for any H. 

(ii) All (up to x <-> y) symmetries (p of the system £^ are 

if = D((p(x, H)) fo r any <p = ((f) 1 , f). 

(iii) In the chosen system of coordinates, the image of the operator □ is closed with 
respect to the commutation in the Lie algebra sym^L- 

(iv) Under a diffeomorphism w = w[w], the r-tuples are transformed by 

Therefore, under any reparametrization u = u[u] of the dependent variables 
u = t (u 1 , . . . ,u m ) in equation £l, and under a simultaneous change w = w[w], 
the operator □ obeys the transformation rule 

u \ w ) w=w[u\ 

u=u[u] 

Consequently, the operator □ becomes well defined: it is a Frobenius operator 
of second kind, see [S|. 

(v) The operator 

A k = n*o (eg)* on (io) 

is Hamiltonian. 
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(vi) The bracket {{ , }}□ on the domain of the operator □ satisfies the equality 

{{>}}□ = {{,}} A k - (11) 

Its right-hand side is calculated explicitly by using the formula ([HI Ej, see 
also [S]) that is valid for Hamiltonian operators = || A"' 3 • D T \\, 

T 

/ r BA a @i \ 

{{p^w Ak = Y,(- i r(D«o [^D T (jpP).-^])(<f<). (12) 

cr,a t,/3 a 

This yields the commutation relations in the Lie algebra sym^L- 

(vii) All coefficients of the operator Ak and of the bracket {{ , }}□ are differential 
functions of the minimal conserved densities w for £l- 

The above theorem is our main instrument that describes all the symmetry genera- 
tors for 2D Toda chains and calculates all the commutation relations in the symmetry 
algebras. 

2. The root system A 2 

Consider the Euler-Lagrange 2D Toda system associated with the simple Lie alge- 
bra s! 3 (C), see (HI HUES], 

£ibda= {u xy = exp(2u-v), v xy = exp(-u + 2v), K=(_^-i)|. (13) 

2.1. The characteristic Lie algebra. First we realize two itegrations of the self- 
adaptive method from [15j, which is based on the use of the characteristic Lie algebra, 
and we obtain two integrals w 1 , w 2 of the characteristic equation D y {w) = on (fl3|) . 
Our reasonings differ from the original approach of [T3] : we do not introduce excessive 
dependent variables and hence do not need to compensate their presence with auxiliary 
constraints. 

Our remote goal is a choice of three layers of the adapted variables 6g, 6q, b{, b\, 
and h\ such that all the coefficients of the linear characteristic equation also become 
linear. Then all the integrals will be found easily, expressed in these variables. The 
number of the adapted variables is specified by the problem, and we have to confess 
that, actually, b\ will be redundant a posteriori because it will be replaced using the 
integral w l in the end. 

Step 1. Regarding the exponential functions 

c(i) := exp ^}' u ' 7 ) 

in the right-hand sides of the Toda equations (CQ) as linear independent, collect the 
coefficients Yi of c{i) in the total derivative 

rn 

D v = J2c{i)-Y i . 

i=l 

Clearly, the solution of the characteristic equation D y (w) = on the Toda chain is 
equivalent to solution of the system {Yi(w) = 0, 1 < i < m}. 
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For system (JT3J), we obtain the vector fields 

d d d 

Y% = h (2u x - v x ) h ( {2u x - v x ) 2 + (2u xx - v xx )) 1 , 

ou x ou xx OU xxx (~\A\ 

3d d 
Y 2 = -x h (2v x - u x ) h ( (2v x - u x ) 2 + (2v xx - u xx )) 1 . 

The underlined terms are quadratic in derivatives of the fields, and it is our task to 
make them linear by introducing a convenient system of local coordinates (see take 2 
of step 3 below). 

Taking the iterated commutators 

Y(ii,...,i k ) - = \Xh >[•••) 5 Yi k ] . . .]] 

of the basic vector fields Y^, we generate the characteristic Lie algebra [TTJ [T5] for the 
Toda chain. If this algebra is finite dimensional (which is the case here), then the ex- 
ponential-nonlinear system (CD) is exactly solvable in quadratures; if the characteristic 
algebra admits a finite dimensional representation, system (CQ) is integrable by the in- 
verse scattering (ibid). For any root system and the Chevalley generators e n , f n , and f) n 
of the semi-simple Lie algebra g, see [5], the characteristic Lie algebra is isomorphic to 
the Lie subalgebra of g generated by the Chevalley generators f n , see [TTj . 
For A 2 , we obtain the commutator 

d d d d 
Y(2,i) = - ~ h t; 3u x h 3^ h ■ ■ 



(This manifests a general fact that is always true: the leading terms of the (k + l)-st 
iterated commutators are the derivations w.r.t. some derivatives u l k+1 , whose order is 
higher than in the leading terms of the preceding, k-th, iterated commutators.) We 
finally note that all the triple commutators, Y(i,2,i) and Y^,2.i), vanish. 

By the Frobenius theorem, a fall of the number of linear independent iterated com- 
mutators at the ith step is equal to the number of first integrals of the characteristic 
equation that appear at this step. The differential order of these new integrals for the 
Toda chains will be i + 1 . 

For the system ([TBI , there appears one (1 = dim(Yj) — dim(F(j l i2 ))) integral, w 1 , of 
order 2. The second and last one (1 = dim(Y(j lii2 )} — (Y^ ui2)i3 ^ = 0)), the integral w 2 , 
has order 3. For arbitrary root systems, the differential orders (shifted by +1) of the 
integrals are described by the proposition in the previous section. 

Step 2. Our remote goal, see above, will be achieved when the expansion 

m m~l 

D x = ^ b i Y i + ^2 & 2 Y (i+i,i) + ■■■ mod 3 : ker D y |^ -»• ker D y \ £ ^ 



is found for the other total derivative, D x . Here the vector field 3 contains only the 
derivations w.r.t. the integrals (yet unknown) and their derivatives, and the dots stand 
for finitely many summands provided that the characteristic algebra is finite dimen- 
sional. The former idea exprimes the replacement of the higher order field derivatives, 
u\ with k ^> 1, using the integrals, while the latter assumption is again based on the 
fact that there are as many integrals as the fields for K Cartan matrices. 
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By definition, put b l := u\. 

Substituting the vector fields Y{ contained in (|14p for d/du % x in D x , we obtain the 
expansion 

D x = u xx Y 1 +v xx Y 2 -\ , 

where the dots stand for the derivations w.r.t. second and higher order derivatives of 
the dependent variables. Consequently, we set 

b\ ■= «L- 

Step 3. Using the four adapted coordinates b l and b{, we rewrite the basic vector 
fields Y k as follows, 

Solving now the system 

Y 1 {w l ) = 0, Y 2 (w l ) = for w\bl b 2 , b{, &?), 
we obtain the integral 

w 1 = u xx + v xx -u 2 x + u x v x - v 2 x . 

Note that, from now on, the coordinate v xx and its descendants can be replaced using 
to 1 , u xx , and first order derivatives. 

Step 1, take 2. Within the second iteration of the algorithm, we repeat steps 1-3 ad- 
vancing one term farther in the expansions. 

Let us indeed replace v xx (although it remains an adapted coordinate) with w 1 . 
Therefore we expand the vector field D y as 

D y = c(l) ■ Y 1 + c(2) ~ mod 3: ker -> 0, 

where the derivations w.r.t. w 1 and its descendants cut off the 'v-part' of the total 
derivative D y . This yields 

9 

UUl XX 

but now the commutator does not contain any derivations w.r.t. the derivatives of v. 

Step 2, take 2. Using the three vector fields, Y\, Y 2 , and Y^ t i), we rewrite 

Dx = u xx Yi+v xx Y 2 +((2u x -v x )u xx -u xxx )-Y {2 ,i)-\ mod3: kerD y \ £ ^ -> ker D y \ ^ . 

Consequently, we set 

b 2 ■ (2u x v x ^)u xx u xxx . 
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Step 3, take 2. Calculating the derivative D y \ £ ^{b\), we substitute it in D y and collect 
the coefficients Yi of the exponential nonlinearities c(i) in this total derivative. 

The result is beyond all hopes: the coefficients of both fields, Y\ and Y 2 , are linear 
in the adapted coordinates, 

In other words, the quadratic terms, which were underlined in (j!4p . are transformed 
into the linear ones. This is due to the quadratic nonlinearity in the new adapted 
variable b\ . 

Finally, we solve the characteristic equation 

Y t (w 2 ) = 0, Y 2 (w 1 ) = for w 2 (bl b 2 , b{, &?, b\) 

under the assumption^ dw 2 /db\ ^ 0. We find the solution 

w 2 = -bl-b 2 bl + blb 2 1 + (blfb 2 -bl(b 2 f. 

Returning to the original notation, we obtain 

+ u x v xx + u x v x - u x v x . 

Obviously, the integral w 2 can be used to replace the derivative u xxx and its differential 
consequences. 

We conclude that now, at the endpoint of the algorithm, both total derivatives, D x 
and D y , contain finitely many terms modulo the vector fields that preserve (respectively, 
annihilate) the kernel ker D y | g . 

In what follows, we do not repeat similar iterative reasonings for the root systems 
B 2 (see ( |T8|) ) and G2 (see p. [HJ), but write down at once the integrals of orders 2, 4 
and 2, 6, respectively. The second integral w 2 for B2 (with a minor misprint in the last 
term) and the higher order 'top' for w 2 for G2 are available in the encyclopaedia [1]. 

2.2. The symmetry algebra: operators and brackets. From the previous section, 
we know the minimal integrals, 

w 1 = u xx + v xx -u 2 x + u x v x - vl, 

UJ ^XXX ^t^x^xx + u x v xx + u x v x - u x v x , 

for the 2D Toda chain (jl~3l) associated with the root system A 2 . Hence we are at 
the starting point for describing its symmetry algebra and revealing the corresponding 
Poisson structures and the KdV-type hierarchies. 
Let us introduce the momenta 

m 1 := 2u T — v T . m 2 := 2v T — u x , 



3 A practically convenient feature of the algorithm is that it allows to fix the 'top' (the higher order 
terms) of the first integrals in advance, whence the redundant freedom in adding derivatives of the 
previously found lower order solutions is eliminated. 
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whence we express the integrals as follows, 



w 1 = 3ml + 3ml - (m 1 ) 2 - t^m 2 - (m 2 ^ 2 



w 2 = 2mL + mL - 2m 1 mi - m 2 m£ + Km 1 ) 3 + |(m 1 ) 2 m 2 - X(m 2 ) 2 - |(m 2 ) 3 . 

By the general scheme of [B], all symmetries (up to x «-> y) of ( 1131) are of the form 
ip = \3((f)(x, [w 1 ], [w 2 ])), where = < (0 1 ,0 2 ) is a pair of arbitrary functions and the 
(2 x 2)-matrix operator in total derivatives is given by formula (Q, 

I — I f^x D x V'x-Dx ~%^x "j^x^x ~%^x ~%^xx 3 V xx\ /I r\ 

U -\ v + £) _l£> 2 +2 w _1 _2 2,2 1 2 I • {^) 

Y ^x \ ^x 3 x ' 3 u xx 3 u xx 3 x ' ^ x x < § u x / 

Next, we calculate the Hamiltonian operator (fTOl) . 

4 - f ^12 
fc \A 21 A 22 

where for the root system A 2 we have that 

A u = 2D 3 X + 2w 1 D x + w l xl 

Al2 = -Dl - w l D\ + (3^ 2 - 2wl) ■ D x + (2w 2 x - wl x ), 
A 21 = D X + w x D 2 x + 3w 2 D x + w 2 x 

A 22 = -\Dl - Iw'Dl - 2w x D 2 x + (2w 2 x - 2w xx - f^ 1 ) 2 ) • D x + \{3w 2 xx - 2w l xxx - 2w l wl). 
The bracket (fT2"j) for A k equals 



rr -. 11 1 1 , T 1 2 211,2/ 22 2 2\ , 2 1/ 2 2 2 2\ 

{{P, 9}}^ = Pz<? - P Qx + [PxxQ ~ P Qxx\ + siP Qxxx ~ PxxxQ ) + 3 W (P Qx - Px<l ) 



(16a) 

{{P. <?}}i fc = [pW - P 1( ll + 2 (pI<1 2 ~ P 2 ^)] + P 2( l 2 xx-P 2 xx ( l 2 - (16b) 

Consequently, not only the image of the entire operator (jToT) is closed under the com- 
mutation, but the image of the first column, of first order, is itself closed under com- 
mutation (see, e.g., [7] or the encyclopaedia pQ, where only the first column of (iT5|) is 
presented). However, the image of the second column of □ = (di,^) is not closed 
under the commutation. Indeed, we box the individual bracket {{ , }}n 1 for the (2 x 1)- 
matrix operator \3i, and we underline the couplings of components in the domain of LI2; 
under commutation, they hit both images of the first and second columns. 

Performing the shift w 2 1— > w 2 + A of the second integral, and taking the velocity of 
the operator Ak, 

4* 



A=0 



we obtain the 'junior' Hamiltonian operator A± = ( 3 ^ 3 ® x ) that is compatible with 
the former. Obviously, the bracket {{ , }} ip) on the domain of A± vanishes identically. 

(We note that the analogous operator = ^| _ (^4fc), where w 1 1— > w 1 + //, is not 
Hamiltonian at all.) 
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The pair (A\,A)-) is the well-known bi-Hamiltonian structure for the Boussinesq 
equation 

w l = 2 ™l - wl x , (17a) 

2 21 211,2 \ 

Indeed, we have that 

5 



to 



| |[^ + |(^) 2 -|(^) 2 + (u; 2 ) 2 



da; = Ak — / w dx. 



Both densities, w 1 and w 2 , are conserved on system ( fTTl) . The symmetry wJ x = 
Afc -p J w 1 dx starts the second sequence of Hamiltonian flows in the Boussinesq hi- 
erarchy 21, see [7] and references therein. 

The modified Boussinesq hierarchy 53 shares the two sequences of Hamiltonians with 
the Boussinesq hierarchy itself, by virtue of the Miura substitution w = w[m] with 
m = m[u}. Namely, for any Hamiltonian ^[w], the flows 

5H[m) SH\m[u]] 

u T = — , m T = L - L 

om ou 

belong to the modified hierarchy <B. The correlation between the two hierarchies, 21 

and 23, and the Hamiltonian structures, 

A k , iS 2) = (A^ 1 , A k , and Bt = (i^)* = B k ,, B k ,, 

for their potential and nonpotential components are standard, see the diagram in [HI 
section 5.1]. The velocities u T constitute the commutative subalgebra of Noether's 
symmetries of the 2D Toda chain ([TBI . 

3. The root system B 2 
The Toda system is specified by the Cartan matrix K = ( _f ~\ ) : 
u xy = exp(2w — 2v), v xy = exp(— u + 2v). 
The integrals for it are of orders 1 and 3 with respect to the momenta: 

w 1 = u xx + 2v xx - 2(v x ) 2 + 2v x u x - u 2 x , (18a) 
w 2 = v Ax + v x (u xxx - 2v xxx ) + u xx v x {v x - 2u x ) (18b) 
+ v xx (Av x u x - 2(v x ) 2 - !//.,)'-) + v xx (u xx - v xx ) 
+ (v x f + (v x ) 2 {u x ) 2 -2{v x fu x . 
Hence the Frobenius operator is 

'u x + 2D X 



v + -D '' 



and 

□i 2 = Dl + (u xx - u 2 x ) ■ D x + {2v 2 x u x - 4:V xx v x + 2u xx v x - 2v x u 2 x + 2v xxx ) : 
□ 22 = Dl + v x D 2 + (2v x u x -v 2 x -u 2 x + v xx + u xx ) ■ D x 

+ (Av 2 x u x - 2v xx v x + 2u xx v x - 2v x u 2 x + 2v xxx - 2v^) . 



KDV-TYPE SYSTEMS FOR RANK TWO SIMPLE LIE ALGEBRAS: AN ILLUSTRATION 13 



The Hamiltonian operator = ( ^j 11 ) has the components 
A u = 10D 3 x + 4:W 1 D x + 2w 1 x , 

A 12 = 3D X + 3w l D 3 x + 6wlD 2 x + (3«4 + 8w 2 ) ■ D x + 6w 2 x , 
A 21 = 3D* + 3w 1 D 3 x + 3wlD 2 x + 8w 2 D x + 2w 2 x} 
A22 = D 7 X + 2w 1 D 5 x + 5wlD x + (Qwl x + Qw 2 + (w 1 ) 2 ) ■ D 3 X 
+ i^ w lxx + ?>w l w l x + 9w 2 x ) ■ D 2 
+ « + 7w 2 xx + (wl) 2 + Aw l w 2 + w l w l xx ) ■ D x 
+ 2- {w l x w 2 + w 2 xxx + w 2 x w 1 ). 

Therefore the components of the brackets ffTTj) for both A^ and □ are 

{{p, q}}h = 2 (p 1( il - pW) + KpIaI - pIqD + (plxQl - pIqD 

+ w1 (pIxqI - pIqIx) + 2 ^ 2 (p 2( i 2 x - pI<i 2 ); 
m q}}h = Hp 2 qI - pIq 2 ) + z{p\ 2 x - pIq 1 ) + z(p 2 qLx - pl xx q 2 ) 

, / 2 2 2 2 \ , r> 1/ 2 2 2 2\ 

+ VPxxQx - PxQxx) + 2w (P Qx - Px9 )• 

Similar to the case of (|TU|) . the commutation of symmetries that belong to the image of 
the second column D 2 of the operator □ = (d 1 , D 2 ) for B 2 hits the image of the first 
column D 1 . 

The 'junior' Hamiltonian operator A\ = -^\ x _ Ak is again obtained by taking the 
shift w 2 t— > w 2 + A in A k : 



AT 



(2) _ f 8D X 



\D X 6D X + 4w 1 D x + 2w} 



The new Hamiltonian operator is compatible with the old one. The bracket on the 
domain of A^ is given by 

{{p, = 2(pW - p 2 q 2 x), {{p, q\} 2 A w = o. 



Likewise to the root system A 2 , the shift w 1 1— ► w 1 + /i produces the operator 
4-\ n Au which is not Hamiltonian. 

The pair (A^\A^\ determines the hierarchy of the KdV-type system 



wl = 6w 2 x , 



w 



D x (2w 2 xx + 2w 1 w 2 ). 



Here we have again that w t = AkT^ J w 2 dx, and the translation w x = At-j^ J w 1 dx 
starts the auxiliary sequence of flows. The construction of the modified hierarchy is 
analogous to the previous case of the root system A 2 , see [H] for details. 
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4. The root system G 2 
The Toda system for K — ( _\ ~\ ) is 

u xy = exp(2w — v), v xy = exp(— 3u + 2v). (19) 

The differential orders of the integrals w.r.t. the momenta equal 1 and 5, respectively: 

w Uxx "i" 3 V xx U x + U X V X 3^ x ; 

w 2 = Uq x - 2u x u 5x + u x v 5x + 10u 4x u x v x - 8u Ax u 2 x - \u ix v 2 x + \u ix v xx 

^V^xUxVx ~\~ ~^V4 X U X gV4 :X V x -\- gV4 x V xx ~\~ ~^V Ax U xx -\- ~^U XXX V XX U X -\- 2,^}u xxx U xx V x 
4:0UxxxU xx U x -\- ~U XXX V XXX ~~^U XXX V XX V X ~\~ ^Qu xxx ll x 1 8u xxx U x V x ~^U XXX V X 

~^U X xxUxV x ~^V XXX U XX V X H ~^V xxx U X xUx 8v xxx U x gV xxx V xx V x ^^xxx^xxUx 

i 26 2 _ q 2,2 3 _i_ J_ 2 __ 17 2 _|_ 40 2 2 _ 98 2 

+ ~3 V xxx U x V x ^ v xxx u x v x ' g V x V xxx i ±8 V xxx ~~3 U xx V xx ' ^U XX U X ZOU xx U x V x 

— 2ll xx V xx v\ + ^§ U xx V x ~ l® u xx ^ u xx v t ~ 5u xx U x V x + 15u xx U x V x — 12lL xx U x V x 
H ^^xx^xxUxVx 3Au xx V xx U x H j^^xx^xx H 3 fii^i^a; 27 ^xx 12l> xx l> x 'U, K 

- f V ^ V l U l ~ l V xx V i + f V lx U l - Av lx U xVx + \v 2 xx vl ~ 3V XX U X - \u x V b x + 

+ §*4^ + fu 2 x v x - 3u 3 x v 3 x . 
The Frobenius operator □ is 

and 



□12 = 2D\ + u x D* + ^15u x u x - \±u 2 x - + 14w xx + 5u xx ) • £> x 

+ (Svxxx ~ 16v xx v x + 10u 2 x v x + ^v xx u x - 44uxxUx - ™u x v 2 x + 24u xx v x - 8u 3 x + 26w a 
+ ^12u^ x - 16v xxx v x - 38u xxx u x - Au 2 x v 2 x + 4m 

xx^ x ^^^xxx^x ^-^-^xx^xx 

+ 2v xx u 2 x + y v xxxU x + Wuxxul + 2lu 4x + 8v 4x - lQv 2 xx - Uu xx u x v x 
+ \v xx u x v x - 51u 2 xx - 9u^j ■ D x 

+ (lSulvl - 23u xx u x v 2 x - y v xxxUxV x + \u x v x + \8u xx u 2 x v x - 12v xx u 2 x v x 
+ ■ 

q ^xx^x^x ^'^^ J xx'^xx^ J x ^-^xx^xx^x 1 ^^xxx^x^x 
+ UxxxVl ~ 20u 4x U x + ^V 2 XX V X + fv 4x U x - 6v 4x V x + 9u 4x V x - 72u XX xUxx 

- 20v X xxVxx + 26v xxx u xx - 8u 2 x vl + 6u xx v x - 8u 2 xx u x - <±v xx vl 



+ 26u XX x v xx "I" 'Z'Uxxx'Ux & u xxx u x 3 V xx U x 9 u x V x 3u xx V x + 14i; xxx Ii x ^ , 
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□ 22 = 3D 5 X + ( 20u x v x - 20u 2 x + fv xx - f v\ + 20u xx ) ■ D 3 r 



+ Uou xxx - 2u x v 2 x - 80u xx u x - y v xxV x + 2u 2 x v x + 4Qv xx u x + ^v xxx + 38u xx v x - \v x J ■ D 
+ {?>8u xxx v x - 25^ - y v xxxV x - 68u xxx u x - 23u 2 x v 2 x + 36ulv x + 9u xx v 2 x 



+ QQu xx v xx + 2v xx u 2 x + 3Av xxx u x + 36u xx u 2 x - 86u 2 xx - 18u x + 34it 4a; 
+ f v ix - AAu xx u x v x + 2v xx u x v x - \v x + 5u x v x J • L> x 

+ ^Gm^m 2 ^ - 22v :Kr:r w. r t>. c - \v\ + 36m^^ + 7u x v x - AQu xx u x v 2 x - 18v xx u 2 x v x 

+ 23v xx u x v 2 x + 40M x:r f - Qu xx v xx v x + Au xxx u x v x + 2Qu bx + f^x 

2 - A0u 4x u x + 6v 2 v x + 20v 4x u x - 9v 4x v x + !8u4 X v x - 120u xxx u x 



^ XXX^ XX I 40'ZJ ^ XXX^ 1 XX 

- 25u 2 x vl + V2u xx v\ - 6v xx v x + 40u xxx v xx + 3v xxx v' x 
,2 



xi 

6 



- 20v xx u x - 18u x v x - Uu xx v x + 20v xxx u x 

Using the fact that the coefficients of the Hamiltonian operator Ak = ^ H ^ 

are differential functions of the integrals u>, we deduce for G2 that! 
A n = fD 3 x + 2w 1 D x + w l x , 

A l2 = 3Dl + fw'Dl + 2 -fwlD x + ((^) 2 + >y • Dl + (f^ - 28wW x ) • D 2 

+ ( 6w 2 + 1|1< - 62w 1 xx w 1 - 62(wl) 2 ) ■ D x + (5w 2 x + fw^ - 32w 1 w 1 xxx - Q&w^wl) 
A 21 = 3Dl + fw'Dl + fw x D x + ((w 1 ) 2 + f O • D 3 X + (9wl xx + 34^X) ■ D 2 X + 6w 2 D x + w 2 x 
A 22 = 2D l x l + 3$w l Dl + \3hwlDl + (414«4 + ^f(w 1 ) 2 ) ■ D 7 X + (819^ + ^w x w x ) ■ D 
+ (1119< + Aw 2 + ^wlw 1 + ±f^) 2 + S^ 1 ) 3 ) • Dl 
+ (10651/4 + + ^f-w 1 w 1 xxx + I22$w l xx wl + eOw^w; 1 ) 2 ) ■ D 4 X 
+ (699w^, + 26w 2 xx + ^wl.w 1 + Sew 1 ^ 2 + ^w l xxx w l x - ^(w 1 ^) 2 - 428w^(w; 1 ) 2 

-88(^.) 2 ^ 1 + 54(^ 1 ) 4 )-^ 3 
+ (303w^ + 29u4 x + ^w\ x w l + 5Aw 2 x w l - b -fw\ x w l x + 54w> 2 - 

- 702w 1 xxx (w 1 ) 2 - WOSw^wlw 1 - 252(^) 3 + 324^(w 1 ) 3 ) ■ D 2 X 

+ (j8wl x + 15w 2 4x + 328w 1 6x w 1 + ^w^w 1 - ^fw\ x w l x + 3Qw 2 x w l x 

- bUwl^w 1 ) 2 + 5Aw 2 wl x - 36w 2 (w 1 ) 2 - 1086(^J 2 - 2000^ 

- 1652«) 2 W 1 - 1120<K) 2 + ^Uw 1 ) 3 + 540K) 2 ( W 1 ) 2 ) ■ D 



5722^1 W 1 

3 ^xxx^xx 



A0 3° W 4:X W xx 



y xxx^x 



4 Here part (|vii[) of our main theorem reveals its true power: a verification for G2 that the bracket 
{{ > }}□> which depends on the fields u through the integrals w l , satisfies equality (JUJ), results in a 
50 Mb size expression. 
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+ 9< + 3w 2 5x + *fv>W + 20w 2 xxx w 



l -fwl x wl + 20w 2 xx wl - ^fw\ xW[ 



- \3&w\ x {w 1 ) 2 + l%w 2 x wl x - l%w 2 x {w 1 ) 2 - §7%w\ x w l xxx - 70%w\ x w l w l x + 18wV 

- 3Qw 2 w 1 x w 1 - w\ xx ■ (U2Awl x w l + 52A(wl) 2 - UA(w 1 ) 3 ) - 680(^J 2 ^ 



Next, we calculate the bracket ffTTj) on the domain of f[2"Uj) using formula (fT2l : 



{{Pi <?}}□ = [PW ~ P l ll] + [f(plx<l 2 ~ P 2( lL) + f(p\x<l 2 x ~ pIQAx) + T iplxxQlx ~ PlxVlxx) 

+ HpLxVL ~ pIxQLx) + ^w l {j> 2 q l xxx - pl xx q 2 ) + 34w 1 (p 2 x q 1 xx - p l xx q 2 x )] 

+ HpIxQ 2 - P 2( lL) + siplxvl - pIqIx) 

j_ 140 1/ 2 2 2 2 \ i of 2 2 2 2 \ i cf 2 2 2 2 \ 

-+- — w (J^y p y 7: J + 6 \P7x1xx PxxItx) i v\P xxx q§ x PgxIxxx) 
+ 376^(p 2 x g 2 - p 2 g 2 J + fw 1 ^ 2 ^ - p\ x q 2 x ) + f u/faL^L ~ pLsL) 

+ ^ '''M/t'/; - pUI) + 2 ^lMxQ 2 - p 2 qD + iss^ 1 ) V<4 - pIq 2 ) 

i 248 1/ 2 2_2 2\, 640 1 f 2 2 _ 2 2 "\ _|_ 4184 1 (2 2 _ 2 2 \ 
' 3 W x\Pax ( 1xx PxxQix) < 3 W xx\PAx < 1x Pxl4x) ' 3 W :ra:a: U^z 9 P ?4xJ 

+ 44w 1 (#L£L - ?4?Lz) + 672w 1 ^(p 2 g 2 :E - p^g 2 ) + l7Q{w l ) 2 {p 2 x q 2 4x - p\ x q 2 x ) 



+ 452(^) 2 (p 2 g 2 ra - pIxxQ 2 ) + ^(w'fipl^q 2 - p 2 q 2 xxx ) + lSu; 2 ^ 2 ^ 2 - p 2 q 2 xxx ) 

+ 26 (™T (pL^L ~ pIxxQL) + 1 ^ 2w1 xx w1 {P 2( 1 2 xxx -PIxxQ 2 ) 

+ 576w 1 ^(p 2 g 2 xa; - p 2 xxx q 2 x ) + ^fw\ x {p 2 xx q 2 x - p 2 x q 2 xx ) + 1360W 1 w\ xx {p 2 q 2 xx - p 2 xx q 2 ) 

+ l ^2wl x wl{p 2 ql x - p 2 xx q 2 ) + 648^(w 1 ) 2 (p 2 ;c g 2 - p 2 q 2 xx ) + 68(^) 2 (p 2 g 2 x - p 2 xx q\ 

+ 772wl x (p 2 xx q 2 - p 2 q 2 xx ) + 36w 2 x (p 2 xx q 2 - p 2 q 2 xx ) + 36(w 1 ) 3 \p 2 xx q 2 x - p 2 x q 2 xx ) 

+ 5Mwl x w\p 2 x q 2 xx - p 2 xx q 2 x ) + 7 -fw\Mq 2 - p 2 q 2 x ) + 1020(«4) Vg 2 - pW) 

+ ?>%w 2 xx {p 2 x q 2 - p 2 q 2 x ) + 1360^4^ (p 2 g 2 - p 2 x q 2 ) + Q%0w\ x w l {p 2 q 2 x - p 2 x q 2 ) 

+ 36^ Vg 2 - p 2 x q 2 ) + 648<> 1 ) 2 (p 2 g 2 - p 2 g 2 ) + 648(^.) 2 ^fe 2 - pVx); 

{{p, <?}}□ = p 2 ^ 1 - p l q 2 x + 5(pIq 2 - p 2 qD + ^(pIxQ 2 - p 2( ilx) + (plxQlxx - pLxqD 

+ 2<W(pLr<? 2 - P 2 qlxx) + ^w l x {p 2 xx q 2 - p 2 q 2 xx ) + 2w\p 2 xx q 2 x - p 2 x q 2 xx ) 
+ 38wUp 2 x q 2 -p 2 q 2 x ) + 18( W 1 ) 2 (/g 2 - p 2 g 2 ). 



hold, where {{p, q}}n = {{&<!}} A. ^ or an y Pi Q( x i M)- This yields all the commutation 
relations between symmetries ip = □(•) of the 2D Toda chain ( fl9l) associated with the 
root system G2. 





The equalities 



and 



[A k (p),A k (q)} = A k (z Mp) {q) - fi^^^) + {{p,q}} Ak ) 
p(p), D(q)] = □(£ D(p - ) (g) - fin®® + {{p, g}} D ) 



(21) 
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Finally, we pass to the KdV-type hierarchy. The deformation ^| A _ Afc under w 2 i— ► 
w 2 + A determines the 'junior' Hamiltonian operator 

m = ( o 6D X 

1 1^6^ AD X + Z^Dl + 54u£D* + (54«4 - 36(w 1 ) 2 ) • D x + 18«4 X - 36u>X, 
It is compatible with the 'senior' operator A^, the bracket on its domain is given through 

{{p, q\Y A m = 36wl {pW - p 2 il) + 18 (AL* - pLxI 2 ) , {{p, q}}%2) = o. 

Applying the Hamiltonian operator At to X J u> 2 dx = (°) , we obtain the KdV-type 
system 

w\ = D x (f < - 32w 1 w 1 xx - 32(^) 2 + 5^ 2 ) , 

w 2 = D x (3wj x + 9wl x + 20w 1 w 2 xx + 18wl x w 2 - ISiw'fw 2 + 

- WwLwl - 7 -fw\ x wl x - tfiwlJ - lSS^l^u; 1 ) 2 - A32wl xx w l x w l 

- 496«) 2 W 1 - 92«4K) 2 + 144<( W 1 ) 3 + lOSK) 2 ^ 1 ) 2 ) . 

(The second equation in this system can be simplified by adding to w 2 a scaling- 
homogeneous differential polynomial in w 1 and thus cancelling some irrelevant terms.) 
The Hamiltonian J w 1 dx starts the auxiliary sequence of flows by the translation 
along x. The corresponding modified KdV-type hierarchy, and the Hamiltonian struc- 
tures B\ = B± , B^i, and By for it, are introduced in a standard way [8]. 

Discussion 

The geometric method for derivation of completely integrable KdV-type hierarchies, 
which is illustrated in this note, is the most straightforward and efficient, to the best 
of our knowledge. Unlike in the fundamental paper |4j, which is based on algebraic 
considerations, the KdV-type systems are derived here, from the very beginning, in the 
bi-Hamiltonian but not in the Lax form. Second, the formalism of pseudodifferential 
operators is not required. We pass from the root systems directly to the 2D Toda 
chains and then to the Poisson structures, skipping over the matrix representations of 
the semi-simple Lie algebras q and the techniques for fixing the gauges. A posteriori, 
the Lax pairs for the KdV-type systems, which are bi-Hamiltonianj w.r.t. the operators 
Af ) and A k , can be obtained in a standard way, Rt + [£f, R] — 0, by using the recursion 

R = Ak o fvii ) and the linearization £p of the right-hand sides F of the KdV- 
type systems w t = F. We argue, however, that this may not be a drawback of our 
concept. Indeed, for q semi-simple, the ambient r-component 2D Toda chain is exactly 
solvable. Therefore, in principle, the modified KdV-type flows should be lifted first 
to the bundles with 2r-dimensional fibres, whose sections f(x) = (/i, . . . , f r ), g(y) = 
(gi, . . . ,g?) determine the general solutions of the 2D Toda chains. We see that the 
liftings of the mKdV-type hierarchies determine the evolution of these Cauchy data. It 
is well known that the Krichever-Novikov equation appears in this context for the root 



Let us finally remark that k may be greater than 2: for instance, we have that k = 3 for the 
Kaup-Boussinesq system, see [5], and two 'junior' structures precede it. 
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system Ai. In our opinion, the posing of the problem of integrability by the inverse 
scattering is more appropriate for that class of evolution equations, which may be not 
bi-Hamiltonian. 

The Lax approach of [I] becomes truely inevitable if the algebras g are Kac-Moody 
algebras, and the Cartan matrices K become degenerate. In this case, our cut-through 
does not work without serious modifications. This will be discussed elsewhere. 

There is one more thing that we lack, it seems, due to the conscious refuse of the 
algebraic language in favour of the geometry of PDE. Namely, this is an explanation of 
the way the 'junior' Hamiltonian operators are obtained through the deformations of 
the canonical operators A^, see ffTUj) . The preference of certain shift directions for the 
integrals must have some explanation in terms of the cohomology of the IV-algebras 
for the KdV-type equations at hand. Also, we do not know a group-theoretic origin for 
the resolvability of the Magri schemes, that is, the existence of the next Hamiltonian 
at each step (equivalently, the vanishing of the first Poisson cohomology with respect 
to the 'junior' operator A\ ). The most prominent technique for proving that is the 
construction of Gardner's deformations for the Drinfel'd-Sokolov systems. 
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